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, Abstract. We analyze in detail three classes of isomondromy deformation problems associ- 

' ated with integrable systems. The first two are related to the scaling invariance of the n x n 

^\ , AKNS hierarchies and the Gel'fand-Dikii hierarchies. The third arises in string theory as the 

representation of the Heisenberg group by [(L'=/")+ , L]= I where L is an n order scalar dif- 
fH , ferential operator. The monodromy data is constructed in each case; the inverse monodromy 

' problem is solved as a Riemann-Hilbert problem; and a simple proof of the Painleve property 

is given for the general case. 
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Introduction 

It has been known since work of Ablowitz and Segur [ASl] that there is an intimate 
relationship between equations like KdV which are integrable by the inverse scattering 
method and Painleve equations; see also [AS2]. It was observed by Flaschka and Newell 
[FN] that just as KdV gives an isospectral flow for the Schrodinger operator, Painleve 
equations are monodromy preserving flows for linear systems with irregular singular points. 
Certain of these problems have been investigated in detail, for 2x2 systems and second 
order scalar problems: [FN], [FZ], [IN], [JMU], [JMl], [JM2]. 
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In this paper we discuss general isomonodromy problems associated to n x n systems and 
higher order equations. The corresponding isomonodromy equations are generally of order 
greater than 2, so they are not the classical Painleve transcendents. However they do have 
the Painleve property, in fact the stronger property that any solution has a single-valued 
meromorphic extension to the entire plane; of. [Ma], [Mi]. We hope, among other things, 
to simplify and clarify the treatment of isomonodromy deformations and their relation to 
Riemann-Hilbert problems, on the principle that the more general the case the less reliance 
on special features. 

The paper is organized as follows. In §1 we describe the formal connection between 
matrix isomonodromy equations and overdetermined n x n systems in two variables. This 
connection is made rigorous in the next three sections, which describe the forward problem 
at the singular points z = oo and z = 0, connect it to a Riemann-Hilbert problem, and 
prove the Painleve property. A few examples of equations and systems which occur in this 
context are: the Painleve II equation 

4:{xu)a; + Uxxx - Qu^Ux = 0; 

the system of three equations of order one 

{xui)a: = aiU + biUjUk, {ij,k} = {1,2,3}; 

the system of two equations of order two with cubic nonlinearity 

{xUi)a; + - u\u2 = = {xU2)x " \u2xx + ului\ 

and the system of two equations of order two with quadratic nonlinearity 

{xUi)x + -^Uixx + '2-U2U2X = = {xU2)x - -^'^2xx + ^UiUix- 

In §5 the rational solutions of the isomonodromy equations are constructed by solving 
finite linear systems. This extends results of Airault [Ai], who found Backlund transforma- 
tions giving rational solutions of some Painleve equations. We develop the gauge theory 
of Backlund transforms in §6. The gauge transformations take the wave functions for one 
solution to those of a new solution and thus transform solutions to solutions. 

Isomonodromy deformations arise from integrable systems in two ways. Some can be 
obtained as self-similar solutions of the given nonlinear evolution equations; see [AS2]. 
This construction is given in §7 for isospectral deformations of n x n first-order operators 
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(AKNS-ZS systems); each of the examples above is of this type. In §7 we also treat the 
isospectral deformations of an n-th order scalar differential operator (the Gel'fand-Dikii 
hierarchy). Examples include the equation 

{xu)x + U+ \Uxxx + \uUx = 

which corresponds to self-similar solutions of the KdV equation and the system 

which corresponds to self-similar solutions of the Boussinesq system. In all these cases the 
Lax pair for the integrable system can be rescaled to obtain a Lax pair for the corresponding 
isomonodromy deformation problem. This was done in detail in the 2x2 case by Flaschka 
and Newell [FN] for the modified KdV equation and its associated isomonodromy problem, 
the Painleve II equation, as well as the sine-Gordon equation and the Painleve III equation 
(cf. also Its and Novokshenov [IN]). In section 8 we show that all self-similar solutions of 
Gel'fand-Dikii flows are in fact solutions of isomonodromy equations, and we treat the 
direct and inverse problems for these equations. The results are analogous to those of §§2, 
3, 4. 

A second class of isomonodromy problems was obtained by M. Douglas [Do] in recent 
two-dimensional theories of quantum gravity. These problems are obtained by replacing 
L by hi in the Lax equation L = [{L'^^'^)-^. , L] for the Gel'fand-Dikii flows. This device 
leads to a representation of the Heisenberg group and to a class of isomonodromy problems 
different from those obtained by the scaling invariance. Two of the simplest examples are 

'^'^xxx ~l~ 2^^x ~l~ ^ — 0, 

j-gti'- -|- gUUxxx ~l" 'g'^xli'xx ~l g"'^^ "^x + ^ = 0. 

This class of isomonodromy problems was treated by Moore [Mol] , [Mo2] . We thank him 
for useful discussions of the topic. In §9 we focus on some aspects of the mathematical 
analysis of these problems and give a different, more complete and self-contained treatment. 

1. Overdetermined systems and isomonodromy equations 

In this section we outline the formal connection between certain overdetermined linear 
n X n systems and monodromy preserving equations. 

Let J and n be diagonal matrices belonging to the space Md{C) of d x d complex 
matrices; we assume that each has distinct diagonal entries and has trace zero. Suppose 
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that q — q{x) is an ofF-diagonal matrix- valued function defined on some real interval or 
some connected complex domain /. Consider the overdetermined system for a matrix- 
valued function il){x^ z) of real or complex x and complex z, i/j & SL{n, C), n > 2: 



(1.1) 
(1.2) 



[zJ + q{x)]'ip; 



dx 

z^ = A{x,z)'ip = ^Z^Aj{x)'lp , An = 11. 



The compatibility condition ( "zero curvature condition" ) for these equations is 

(1.3) 
Set 



[^-(zJ + q),z^^-Aix,z)]=0. 



(1.4) 



Fn-j, n> 1; Ai = F^-i + xJ; Aq = + xq{x). 



Then the zero-curvature condition (1.3) is equivalent to the sequence of conditions 



(1.5) Fo = 

together with the equation 
(1.6) 



d 



[ ^ - ? , a;? F„ ] =0 



The conditions (1.4) do not determine the Fj uniquely but, as we shall indicate, there 
is a determination for which the Fj are polynomials in q and its derivatives. Then (1.6) is 
a nonlinear differential equation for q, and we show that it is an isomonodromy equation 
for the system of differential equations in the complex domain (1.2): under the flow (1.6), 
the monodromy of (1.2) is invariant. 

Examples: 1. The simplest nontrivial example is obtained with n = 2, k = 3, ij, = J = 
diag(l, —1), and 

q{x) = u{x) 



1 

1 



Then 



Fi = q, F2 





—V? Ux 








^ 2 


—Ux 


' 4 



— 2u^ 

Ur.r — 
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and equation (1.6) becomes the equation satisfied by self-similar solutions of the mKdV 
equation, namely the Painleve II equation 



See §§7, 8 for a full discussion of self-similarity. 

2. Although we assume generally that n is at least two, certain cases with n = 1 can be 
treated, e.g when J and n are both real. If we assume this and assume that q + = 0, 
then the 3x3 case with n = 1 and k = 1 leads to a system of equations of the form 



where the a^'s and 6j's are constants. These equations correspond to self-similar solutions 
of the three-wave interaction equation. Some more examples are sketched at the end of 
this section. 

We suppose from now on that q is a smooth function on the domain / C R and that it 
is off'-diagonal: 



4:(xu)x + Uxxx - Qu^Ux = 0. 



{xui)x + aiu + biUjUk = 0, 



{z,j,/c} = {1,2,3}, 



q{x)jj = 0, 



l<j<d. 



Let 



(1.7) 



il>{x, z) = m{x, z)e 



^{x, z) = -Z^jJL -\- XZJ. 



The equations (1.1), (1.2) for are equivalent to 



(1.8) 



(1.9) 



dm 

dx 
dm 



A{x^ z)m — m {z^n + xzJ) . 



[zJ , m] + qm , 




By a formal solution of (1.8) we mean a formal power series in z 



-1 



oo 



m 



{x,z)=Y,z-^f,{x), fjeC^{I;Md{C)), fo = l 



which satisfies (1.8) formally; this means that for each N 
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Equivalently, 

(1-10) £-^-^^ = ['^'-^^+^]- 

These relations determine the off-diagonal part of fj+i from fj by inverting ad J. At the 
next step they determine the diagonal part of fj+i, up to a constant matrix, in terms of 
the off-diagonal part of fj+i since the diagonal part of [ J , fj+2 ] vanishes and the diagonal 
part of qfj+i involves only the off-diagonal part of fj+i- Thus these relations are solvable 
recursively, so formal solutions exist. Note in particular that q = fiJ — Jfi and that the 
diagonal part of /i may be taken to be 0. For later use we need to go one step further 
in this discussion. Let us ftx a point xq and choose the unique formal solution with the 
property that fj{xo) is off-diagonal for all j > 0. Then one can prove by induction that 
each Taylor coefficient at xq of each fj is given by a universal polynomial in the Taylor 
coefficients at xq of the entries of q. Moreover, each of these polynomials has constant 
term 0. In particular, fj+i{xo) is a polynomial in the j-jet of q at xq. 
We define formal solutions of (1.9) in an analogous manner. 

Theorem 1.1. Let m — Sj^o "'/i(^) ^ formal solution of (1.8). The formal series 

CO 

(1.11) F = mum-^ = ^^z-^Fjix) 

j=0 

is independent of the choice of m. The coefficients Fj = Fj^f^ are traceless polynomials in 
q and its derivatives: 

(1.12) Fj{x) = PM^), Q'i^), Q'~\x)). 

Moreover, the Fj satisfy the conditions (1.5). 

If m is also a formal asymptotic solution of (1.2), then the relation between the coeffi- 
cients Aj and Fj is given by (1.4) and q satisfies the isomonodromy equation (1.6). 

Proof. Suppose that mi and m2 are two formal asymptotic solutions; then m2 has a formal 
inverse and g = m]^^m2 is a formal asymptotic solution of the equation dg/dx = [zJ,g]. 
The corresponding relations are 

and since (/o = 1 it follows recursively that each gj is diagonal and constant. Therefore m2 
is a formal power series product m2 = mig with g diagonal and so m2//m2 ^ = mi//m^^. 
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It is clear from the preceding that mjji'm''^ is independent of the choice of formal solution 
m. The assertion about the coefficients of mjj,m~^ may be derived by localizing and using 
a result of [Sa], but we give here another proof which adapts more readily for use in §8. 
Fix any xq in the domain of q, and let m be the formal solution whose diagonal part at xq 
is 1. From the remarks above we conclude that the coefficients of mfim~^ at xq are given 
by universal polynomials in the entries of q and their derivatives at a^o- But this fact is 
independent of the choice of m and of xq. 

If m is a formal solution of (1.8), then F = mjimT^ is readily seen to be a formal 
solution of 

(1.13) [^_^j_,,F]=0. 

Therefore its coefficients satisfy the recursion relations corresponding to (1.13), which are 
(1.5). 

Finally, suppose that m is also a formal solution of (1.9). Then zrug has no constant 
term, so we conclude that (as formal power series in z~^) 

A{x, z) = m{x, z) [z^ + xzJ] m{x, z)~^ + 0{z~^) 

= z'^F + xzJ + xfiJ - xJfi + 0(z-'^) = z'^F + xzJ + xq + 0{z-^) . 

This is equivalent to (1.4). Now (1.3) can also be obtained by equating terms in formal 
power series expansions, and we deduce the isomonodromy equation (1.6). □ 

Let us note explicitly that 

Fo = ii, Fi = (ad J)"^[//, g] , 
Fj+i = I — j (ad J) ■'[fJ' , q] + { terms of order < j}. 

We assume from now on that the Fj are those in Theorem 1.2. 

Corollary 1.2. The isomonodromy equation (1-6) is an algebraic differential equation for 
the matrix function q. 

Further examples: 3. Here we take n = 2, J = /j, = diag(l, —1), and k = 2, with q a 
general off-diagonal matrix, 

^0 ui 

U2 
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Some computation shows that Fi = q, F2 — \qJq + ^Jqx, and equation (1.6) becomes 
{xq)x + ^JQxx - Q^Jq = 0. Explicitly, 

{xui)x + \uixx - u\u2 = = {xU2)x " \u2xx + u\ui. 

If we replace J by iJ the formula is the same and is compatible with the reductions 
U2 = iwi, in which case the solutions are exactly the self-similar solutions of the cubic 
nonlinear Schrodinger equations. 

4. Finally we take n = 3, J = diag(Q:, a^, 1), where a = e^'^*/^ is a primitive cube root 
of 1, = J^, and k = 2. We assume that q has the form q = uiU. + U2li^, where ui and 
U2 are scalar and 11 is the permutation matrix 



n 



1 

1 

1 



Then Fi = [J^,r] and F2 = Jtx + TxJ + q — 2uiU2l, where 



r = (1 - a)~^uiJ^U + (1 - a^)~'u2JU 



2\-l„ 



Equation (1.6) becomes {xq)x + Jtxx +fxx + {.(f')x — 2(t(.it(.2)xl, or explicitly the coupled 
quadratic nonlinear system 



{xUi)x + -^Uixx + 2u2U2x = = {xU2)x - -^'^2xx + 2uiUix- 



2. The forward monodromy problem at 2; = 00 

In this section we continue to assume that g is a smooth, off-diagonal matrix function 
on a real interval or a connected complex domain / and we assume it satisfies the algebraic 
differential equation (1.6). (Consequently q is analytic in /.) 

Using a partial sum of a formal solution of the x-cquation as a gauge transformation 
(or parametrix) we show that the z-equation has actual solutions which are well-behaved 
in specified sectors and have a special form. Moreover, the classical Stokes matrices, which 
link the solutions from different sectors, are independent of x. (It is classical to use a 
formal solution of the z-equation as a gauge transformation to obtain a tractable integral 
equation for the transformed solution; cf. Chapter 5 of [CL]. The special feature here 
is that the phase function which occurs in the ^-asymptotics of the solutions has a very 
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special form, and this is a consequence of the fact that the formal solutions are also formal 
solutions of the a;-equation.) 

Let the polynomial part of the formal series z'^F be denoted 

n 

We assume that q satisfies the isomonodromy equation (1.6): 

+ =0. 

In view of the recursion relations for the Fj, an alternative form is 

{xq)a, + [J, Fn+l] =0. 

It is convenient to introduce a condensed notation for the differential operators which recur 
throughout: 

D^ = — -{zJ + q); D, = z—- [z''F]+ - x{zJ + q) 

Then the conditions (1.5) are satisfied and (1.6) is the remaining requirement for the 
commutator condition 

(2.1) [D^,D,]=0. 

We shall show that (1.2) has certain normalized solutions with asymptotic expansions 
in sectors of the plane. To describe these sectors, we set 

T.{3,k) = {zeC■.Rez''{^ij- ^ik) = Q}, l<j<k<d; E=\jE{j,k). 

Each E(j, k) is a union of n lines through the origin, with angle n/n between successive 
rays. The distinguished open sectors {Oj^} which we consider are described as follows. 
Each sector is bounded by two rays from E, and its interior contains exactly one ray from 
each E(j, /c). We take the leading edge of a sector to be the boundary ray with greater 
argument and we number sectors in order of increasing argument of leading edge, from 
some starting point. 
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Theorem 2.1. Suppose that q satisfies the isomonodromy equation (1.6). Then for each 
X and each sector Q.i, there is a unique solution of 

(2.2) = ([z"F]+ + a;zJ + a;g)V', z^n^ 

oz 

which has the form ipi, = mi^e* as in (1.7), where the rui, have the (same) asymptotic 
expansion 

oo 

(2.3) ruj^ ~ '^^z~^ fj{x) as z ^ oo, 
with /o = 1. 

The expansion (2.3) is valid uniformly as z tends to infinity in any closed subsector of 
Qi, U {0} and as x varies in any compact set. 

There are constant matrices (Stokes matrices) Sjj such that the relation between succes- 
sive solutions is 

(2.4) ipyj^-i^x^z) ^ iii,{x,z)Sy, ze9.v^^v+i- 
Finally, the ipy satisfy 

(2.5) det = 1 ; ^ = {zJ + 9) v.. 

The rest of this section is devoted to the proof of Theorem 2.1. We begin with a gauge 
transformation to convert (2.2) to a more tractable form for large z. Let m — fj be 

a formal solution of (1.8) with fo — 1 and set 

n 

f{x,z) = '^z~^ fj{x). 

3=0 

Then fiif-^ = F + 0(^-^-1), while (1.8) implies that {zJ + q)f = zfJ + 0{z-^), so 

(2.6) {[z''F]+ + xzJ + xq)f = f (^"// + xzJ) + 0{z-^) . 

If we look for a solution of (2.2) in the form ip = ftp, then (2.2) is equivalent to 

{z^~^IJ, + xJ) ip + r(a;, z) ip 



dip _ , 



dz 

with r{x, z) = 0{z~'^). Setting ip = rhe^, we convert this to the integral equation 
(2.7) m = 1 + r e*(^)-*(^)r(C)m(C)e-^(^)+^(^) dC, 

J oo 

where we have suppressed the dependence on x. Note that (2.7) is a matrix equation and 
the path of integration may differ entry by entry. 
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Lemma 2.2. In each sector Q^,, the paths of integration in (2.7) may be chosen so that 
there is a solution for large z in each closed suhsector, with an asymptotic expansion valid 
uniformly in the smaller sector. 

Proof. Let 

^xzU, /c) = { C e C : Re[$(a;, z) - ^x, C)]jj = Re^x, z) - ^x, C)]kk }, 3 < k. 

This set consists of n regular curves, one of which goes through z. For z not in E = IJ k) 
the curve through z is asymptotic to the nearest two rays of Tj{i,k). It foUows that for 
each 2 G O;^ we may choose the path of integration for the (j, k) and (fc, j) entries in (2.7) 
to he along this branch and to be asymptotic to the ray of E(j, k) which lies in the interior 
of VLjy. The effect is that conjugation by e**^^)"*''^'' multiplies each entry by a hounded 
exponential. 

If z lies in a closed subsector of 0,^, the associated paths of integration lie in {|C| > R} 
when \z\ is large. Since r = 0{z~'^), it follows that eventually the integral equation (2.7) 
with our choice of contour has a unique solution obtained by successive approximations. By 
a standard argument, this solution has an asymptotic expansion which is valid uniformly 
in each closed subsector. □ 

Note that the solutions so obtained are holomorphic, and thus extend to the entire 
sector; indeed they extend as functions on C \ (multi- valued in general). 

Lemma 2.3. Suppose is a second solution of (2.2) in Q^, for fixed x, with the property 
that the limit of '4^ye~^ as z tends to infinity along any ray in Sy isl. Then V'ly = '4^v 

Proof. There is a constant matrix S such that -01^ = '^i^S'. From the asymptotic expansions 
it follows that the conjugation 

converges to 1 as 2 tends to oo along any ray in 0,^. Then the diagonal entries of S are 1. 
By boundedness, 

= if Re2"(//j - //fc) > 
along such a ray. By assumption, 0,,^ contains such rays for each j ^ k, so S = 1. □ 
Lemma 2.4. The functions ip^, satisfy (2.5). 

Proof. The commutativity property (2.1) implies that the function 

d 
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is also a solution of (2.2), so D^il^i, = ifj^T^x) for some matrix- valued function T. The 
construction of V'lv shows that the corresponding m^y is differentiable with respect to x 
and the asymptotic expansion of rrii, can be differentiated term by term with respect to x. 
Then 

(2.8) m-^(^^ -zadJ-?^m^ = e*Te-*. 

The left side is bounded as 2; — > oo so, as in the proof of Lemma 2.3, it follows that the 
matrix T must be diagonal. Thus the right side of (2.8) is independent of z. On the other 
hand, a check of the asymptotics of the left side of (2.8) shows that the leading term for 
large z is [mi^i, J] — where ttt-i^i is the coeflBcient of z~^ in the asymptotic expansion 
of m^. This term is off-diagonal since q is off-diagonal. (This is the first time we have 
used the assumption that q is off-diagonal.) Thus both sides of (2.8) must vanish, hence 
Dxi^v = 0. We have proved in the process that 

(2.9) q{x) = —[J ^ rrii^i]^ where rrii, = 1 + z~^mi,i + 0{z~'^). 

Finally, the equation (1.8) for m,^ implies that detm,^ is constant, hence is identically 1. 
We have assumed that $ is traceless, so ipy — rriye^ also has determinant 1. □ 

Lemma 2.5. There are constant matrices Si, such that for all x in I 

(2.10) ipu+i{x,z) = 'ipu{x,z)S„, 2; e fiix+i n n^y. 

Proof. This is clear from the fact that in the common domain of definition, both functions 
satisfy linear first order systems of differential equations in x and in z. 

The preceding lemma makes possible the final step in the proof of Theorem 2.1. 

Corollary 2.6. The coefficients of the asymptotic series for ruy do not depend on the 
sector. 

Proof. By Lemma 2.3, 

$ o — 1 — $ 

rrii, = m^+ie e 

on any ray in 0,y_|_i flOj^, and our arguments show that the only non-zero entries of S^, must 
remain bounded under the conjugation by as z tends to infinity on the ray. Therefore 
the off-diagonal entries decay exponentially as 2; ^ 00. This implies that the asymptotic 
expansions of m^y and mj^+i are the same. 
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3. The forward monodromy problem at 2; = 

To complete the analysis of the 2;-equation one must examine the behavior at the regular 
singular point z = 0. 

We continue to assume that q is a solution of the isomonodromy equation (1.6) on a 
domain /. We take the sectors and the matrix functions V'l^, m^, $ as in §2, and we 
assume that there are N sectors numbered cyclically: Qn = ^^o- 

Theorem 3.1. There are invertible constant matrices Ci, and a matrix function U{z), 
holomorphic and invertible in C\0, such that 

is regular at z = for every x E I and is independent of v . 

Lemma 3.2. There is a fundamental solution of Dxi/j = 0, Dzi/j = which has the form 

(3.1) ip{x,z) = f{x,z)V{z)z^, zeC\0, 
with f{x,-) entire, V entire, B constant, and detiVz^) = 1. 

Proof. Choose xq E I and let /(x, z) be the unique solution of Dj;f = which satisfies 
the initial condition f{xo,z) = 1; it is entire as a function of z. Since [Dx-,Dz] — 0, 
it follows that D^D^f = D^D^f — 0, so Dzf is also in the kernel of Dx- Therefore 
Dzf = f{x, z)C{z) for some (entire) function C. Note that det / is constant with respect 
to X, hence is identically 1. We look for ip in the form '^(a;, z) = f{x, z)U{z). The necessary 
and sufficient condition for Dz{fU) = is 

(3.2) z'^ = -CU. 

dz 

From the equation D^f = fC and the condition /(xq, ,2) = 1 we find that C is a polynomial 
in z and — C(0) = x^q^xo) + Fn{xo). Equation (3.2) has a regular singular point at the 
origin, so it has a fundamental solution U {z) = V{z)z^ , where V is entire and B is constant 
[CL, Ch. 4, Theorem 4.2]. Now det U is constant, so we may choose det U = 1. □ 

Proof of Theorem 3.1. Each fundamental solution ip^, = m^e^ differs from the solution 
constructed in Lemma 3.2 by right multiplication by an invertible constant matrix Cy. □ 

Generically one has more information about the fundamental solution at 2 = than is 
given above. We begin with a general remark about the constant term Aq{x) — xq{x) + 

Fn{x). 
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Lemma 3.3. The isomonodromy equation (1.6) implies that the matrices xq + are 
similar for all values of x. 

Proof. Let a{x) be a non-singular solution of da/dx = qa. Then 



d 

dx 



a ^{xq + Fn)a 



= a 



d 

q,xq + Fn 

dx 



a = 



so a ^{xq + Fn)a is constant. □ 

Corollary 3.4. Ifq and its derivatives of order less than n have limit zero anywhere, then 
xq + Fn vanishes identically. 

Definition. The matrix function q is proper if -D^V' — has a fundamental solution of the 
form 



(3.3) 



ip{xQ^ z) — w{xq^ z)z'^°^^°\ w{xo,-) entire, w{xo,0) 



for some xq G /. 

Remarks. This is a condition on the (n — l)-jet of q at xq. By the proof of Lemma 3.2 
the condition is in fact independent of xq, and applies to all x G /. 

In particular, q is proper if for some xq no two eigenvalues of A{xo) differ by a non-zero 
integer (CL, Theorem 4.4.1]); hence generically, q is proper. 

When the eigenvalues of Ao{xo) differ by an integer, there is still a solution at the 
origin of the form wz^, with w entire ( [CL], Ch. 4, Theorem 4.2); but generically, w 
is degenerate at the origin, and B need no longer be similar to Ao{xo). {B is similar to 
Aq{xq) if w is invertible at the origin.) 

These two situations are simply illustrated by the Painleve II equation. 

Example: Half-Integer Solutions of Painleve II. As in §1, we take 



Fi = q, 



J = II 



1 
-1 

-u^. u'^ 



ua = u 



Ua 



1 

1 

Uxx - 

-2^3 



As noted earlier the isomonodromy equation is 



(PII) 



\Uxx - +XU^U 
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where f is a, constant; and Ao{x) = va is constant. 

Let us consider the case v =1/2 and look for a solution of the form wz^l"^ . We pick an 
initial point and let 

u{xq) = a , u {xq) = (3 . 



The equation Dz{wz'^^'^) = is 



z^iwz"''^) = zw^z''!'^ + wz'^l'^o = (Ao + + . . . )wz''l'\ 
dz 

Substituting a power series expansion w = wq -\- wiz -\- . . . into this equation, we get the 
recursion relations 



/2. 



r/2 



where 



Ai = F2 + xq = 



xo - o?l2 13/2 
-/3/2 -xo + a^/2 



From the first equation, wq commutes with u; and since the equations are invariant under 
right multiplication by a constant matrix, we may factor out wq on the right. Hence 
without loss of generality we may take wq = 1. Taking 



Wi 



(I b 

c d 



and turning to the second recursion relations, we find a + d = b + c = and the constraint 



(C) 



P + a'^ = 2xo. 



Since this constraint is independent of the initial point xq we get the additional equation 



u' + u'^ = 2x. 



It is easily seen that any solution of this equation is also a solution of (PH). Moreover, 
setting u = xVx one finds that x satisfies the Airy equation 

X" - 2xx = 0. 

These special solutions were originally discovered by Gambler [Ga]; other special half- 
integer solutions were constructed by Airault [Ai], using Backlund transformations. 
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The constraint C on the 1-jet of the solution at xq shows that the proper solutions of 
Painleve II for v =1/2 form a submanifold of codimension 1, and so are non-generic. The 
same considerations apply in general to solutions of the isomonodromy equations when the 
eigenvalues of differ by integers. The monodromy data for the general case is discussed 
in Theorem 4.1. 

It is important to note that q need not be proper. The condition that it be proper 
is a non-trivial condition on the (n — l)-jet of g, which determines the operator D^; see 
Chapter 4 of [CL]. 

Theorem 3.5. Suppose that q is proper. Then there are unique constant matrices B^, 
such that the functions z)z~^'' are entire functions of z, invertihle at z = 0. Each 

By is similar to Aq{xq) and therefore has trace 0. Moreover 

(3.4) = S'^^BySy ; exp(27rii?j,) = S^S^^i ■ ■ ■ S^-i. 

Proof. Combining the assumption that q is proper with the construction in Lemma 3.3, 
we obtain a fundamental solution 

with w{x^ .) entire and t(;(a;o,0) = 1. Then there are constant matrices Ci, such that 

(3.5) = = tua-z^" , By = C-^AQ{xQ)Cy. 

The first relation in (3.4) follows from (2.4). The second relation follows from the fact 
that the analytic continuations of '^i, and z^" around the origin are '^vSvSn^i • ■ -81,-1 
and 2;^" exp{2TTiB,y) respectively. 

If wz^ = wz^ with w and w regular and invertible at the origin, then z^z~^ is regular 
at the origin and it follows that B = B. Thus the B,^ are unique. 
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4. The inverse problem and the Painleve property 

The inverse problem, to determine the solution of the isomonodromy equation from 
its monodromy data, can be formulated as a Riemann-Hilbert problem. This idea was 
introduced in the context of evolution equations by Shabat and utilized by Flaschka and 
Newell in 2 X 2 cases. We give the general formulation in this section. We also show that 
every solution of a Riemann-Hilbert problem of this type is associated to a solution of an 
isomonodromy equation (1.6). 

It has been shown by Malgrange [Ma] and Miwa [Mi] that equations like those considered 
here have the Painleve property. We sketch here an argument based on the Riemann- 
Hilbert problem. 

Suppose again that qis solution of the isomonodromy equation. 

Theorem 4.1. The function q is uniquely determined by its Stokes matrices Sv, the func- 
tion ^{x^ z) = ^-z"a* + xzJ, and the functions Uj^ of Theorem 3.1. 

If q is proper, then q is uniquely determined by the Stokes matrices Si,, the function ^, 
and the exponents By of Theorem 3.5. 

Proof. Let the functions 1})^, rrij^ be as defined previously. Let be a ray in the intersection 
Vty n Oj^-i-i, let r = = 1} be the unit circle, and let 

K = V\J (|JS^n{^ : \z\ > 1}) . 

Fix a; e / and define a function M on C \ K by 

(4.1) M{z) = i)y{x, z)Uy{z)-^ , for \z\ < 1 

(4.2) M{z) = my{x,z) for 1^1 > 1, z lying between "E^-i and . 

(Recall that by Theorem 3.1 the value given by (4.1) is independent of u.) The matrix 
function M has the properties: 

(4.3) M is holomorphic and invertible on C\K ; 

(4.4) M has limit 1 as 2 — > 00 ; 

(4.5) M is continuous up to the boundary from each component . 

Let Mr and M^, denote the boundary values of M on the circle from the unit disc and on 
the boundary of the region outside the circle between the rays and T,,^, respectively. 

These boundary values are linked on their common domains of definition by 

(4.6) My{z) = Mr(;2)C/,e-*(^'^) ; M,+i(^) = M^e^^^'^^^^e"*^^'^) . 
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Uniqueness in the general case will follow if we show that M is uniquely determined by 
the Riemann-Hilbert problem (4.3)-(4.6). But this is standard: if M is a second solution 
of (4.3)-(4.6), then MM~^ is piecewise holomorphic and continuous, with value 1 at oo, 
so MM-i = 1. 

Suppose now that q is proper. Then we redefine M on the unit disc by 

(4.7) M{z) = ip^{x,z)z-^'' if 1^1 < 1 . 
Then M has the previous properties, with (4.6) modified to 

(4.8) M^{z) = Mr{z)z^^e-'^^''''^ ; M^+i{z) = M^e^^^'^^^^e'^^^'^) . 

As before, this M is uniquely determined by the conditions (4.3)-(4.5) and the boundary 
relations (4.8). □ 

Our principal result concerning the isomonodromy equation is the following. 

Theorem 4.2. Any solution of the isomonodromy equation (1-6) on an interval or con- 
nected domain has a single-valued meromorphic extension to all x e C. In particular (1-6) 
has the Painleve property: the only movable singularities of its solutions are poles. 

Proof. The Riemann-Hilbert problem (4.3) to (4.6) has data which depend holomorphically 
on X through the function $(a;, z) — $(0, z) + xzJ. 

Since the ray E,^ lies in the sector fit/, the off-diagonal elements of e^S'j^e"* decay 
at infinity like exp(— el^l"^), which more than offsets any growth from exp(a;^J) and its 
derivatives with respect to x. 

The modified Riemann-Hilbert problem can be written as a singular integral equation 
with parameter x, which we write symbolically as 

(4.9) M = 1 + M , 

where Cx also depends on the data 5", U . The operator Id — Cx depends holomorphically 
on x in the entire plane and is invertiblc for all x in the original domain / of q. It is 
shown in [BC] that Id — Cx is Fredholm with index zero; in fact it is shown that inverting 
Id — Cx can be accomplished (if at all) in two steps: inverting a small perturbation of 
the identity, followed by solving a finite system of linear equations which depend on x. In 
the present case both the perturbation of the identity and the system of linear equations 
depend holomorphically on x, and it follows that the solution is meromorphic with respect 
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to X in the entire x-plane. Therefore the associated potential q is also meromorphic in the 
entire a;-plane. (In [BC] the sufficient conditions include compatibility conditions for the 
data at the intersections points of its domain of definition; these conditions are shown in 
[BDT] to be a consequence of a product condition on Taylor expansions of the data at the 
intersection points. In the present case the data is assumed to have come from the direct 
problem, so the conditions are automatically satisfied.) □ 

The next result closes the circle, insofar as proper solutions are concerned. Recall the 
following facts about the Stokes matrices and the matrices Bi,: 

(4.10) e'^S^e~'^ is bounded as 2 ^ cxo, z G , {Su)jj — 1 ; 

(4.11) tr(S^) = exp(27riS^) = S^.S^.+l ■ ■ ■ S^-i . 

Theorem 4.3. Suppose matrices Si,, Bn satisfy the conditions (4.-10), (4-11)- If the 
associated Riemann-Hilhert problem (4-3)-(4-5), (4-8) has a solution M{xq,z) for some 
value xq of X, then it has a solution M{x,z) except for a discrete set of values of x & C 
For \z\ > 1, z ^ [J Til,, the function ip{x, z) = M{x, z)e^ satisfies the system of equations 

(4.12) D^ip = 0, D,ip = 

where the operators = d/dx — zJ — q and = zd/dz — [z'^F]^ — xzJ — xq are defined 
by a unique off- diagonal matrix function q. The function q satisfies the isomonodromy 
equation and is proper. 

Proof. Because of the properties of the Riemann-Hilbert problem as indicated in the pre- 
vious proof, solvability ai x = xq implies solvability except in a discrete subset of C. The 
determinant satisfies a Riemann-Hilbert problem which forces det M = 1; in fact since 
tr($) = tT{B,j) ~ 0, it follows that detM is continuous across K and hence entire in z, 
and (4.4) implies that det M is 1 at infinity. Thus M is invertible. The function 

^ -^adJ ) M 
ox J 

is piecewise holomorphic, continuous, and bounded, hence has value q{x) independent of 
z. Therefore i/j = Me* satisfies the first of the equations (4.12), for this choice of q. 
Moreover M{x, ■) has an asymptotic expansion J2 fji^) as 2; — > 00 and q = —[J, fi] is 
off-diagonal. On the other hand, the function 



d 

z— z'^ aAu — xzaAJ ] M 

oz 



M 



-1 
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is piecewise holomorphic, continuous, and 0{z'^ ^) as 2; — > 00, hence is a polynomial of 
degree less than n in z. Therefore ip satisfies an equation 



As in the proof of Theorem 1.1, we conclude from the asymptotic expansion of M at 
infinity that 



Thus ip satisfies the equations D^ip = and Dgip = 0, so q satisfies the isomonodromy 
equation. 

Finally, we must show that q is proper. Let ipy be the restriction of i/j to Q^H {\z\ > 1} 
and extend i/Ji^ to the unit disc by setting ifj^, = Mz^'^ for \z\ < 1. Then i/j^, is continuous 
across the circular arc, because of (4.8), and is therefore holomorphic. Since M is regular 
at the origin, ipi, has the desired form and q is proper. □ 

5. Rational solutions 

In this section we investigate rational solutions of the isomonodromy equations. These 
are the analogues of the refiectionless potentials in inverse scattering theory; indeed there 
is a close connection between the Stokes matrices Si, being trivial (=1), the normalized 
wave functions rui, being rational in z, and q being rational in x. 

Theorem 5.1. Suppose that q satisfies the isomonodromy equation (1.6). The associated 
system (1.8), (1.9) has a solution m which is rational in z and equals 1 at z = 00 if and 
only if each Stokes matrix Si, equals 1 . 

// also q is proper, then q and m are rational functions of x. 

Proof. If (1.8), (1.9) has a solution m which is rational in z and equals 1 at 00, then by 
uniqueness rrii, = m, all and so each 5",^ = 1. Conversely, suppose each Si, = 1. Then 
the nil, coincide and so m = nii, is regular on the Riemann sphere minus the origin. It 
follows that in the representation (3.1), m = /(x, z)V{z)z^ with fV entire, the factor z^ 
is single- valued. Therefore exp(27rzi?) = 1, and therefore B is diagonalizable, with integer 
eigenvalues. This shows that the origin is a pole for m, hence m is rational. 

Now suppose that q is proper, i.e. that B has trace zero. Suppose that m has a pole 
of order r at ^ = 0. Then the smallest eigenvalue of B is — r. Denote the dimension of 
the eigenspace corresponding to the eigenvalue s by ds- Let {vi,V2, ■ ■ ■ , Vd) be a basis of 




A{x, z) = + xzJ + xq . 
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eigenvectors of B with eigenvalues s = ai > a2 > ■ ■ ■ > cxa = —r. For 1 — r<k<s = ai 
define a d x (dk + ■ ■ ■ ds) matrix Vk by taking the columns to be those vj with aj > k. Let 

oo 

where the matrix-valued function Cj is a polynomial of degree j in x. Then the condition 

(5.1) me^z-^ = [l + z-'^fi + --- + z-''fr]e^z-^ is regular at ^ = 
is equivalent to the sequence of equations 

(5.2) /,Vi_, = 0, 

[/rCi + fr-l] V2-r = 0, 

[f^Cr-l + --- + fl]Vo = Q, 
[frCr + ---+flCi]Vi^-Vi, 

[frC,+r-i + --- + hCs]Vs = -Vs. 

The first equation in (5.2) represents d{di-r + d2-r + - ■ ■ ) linear equations for the d'^ entries 
of fr- Altogether the number of linear equations represented by (5.2) is 

d [{di-r+d2-r + ...) + (rf2-r + C^3-r + •.•) + (c^S-r + ...) + ...] 
= d [di-r + 2d2-r + 3(i3-r + . . . ] 

= d[{l- r)di-r + (2 - r)d2-r + ...] + dr [di-r + c^2-i + • • • ] 

= d [tr(S) + rd-r] + dr[d- d^r] = rd^- 

Thus (5.2) is a set of rd'^ equations for the rd'^ entries of the fj. As noted above, the 
coefficients in these equations are polynomials in x. To prove m is rational in x we only 
need to show that the system has a unique solution. By assumption it has one solution. 
As noted, (5.2) is equivalent to (5.1). Note also that 

lim det(/e*2--^) = 1 , 

Z—^00 

SO (5.1) implies that fe^z~^ is invertible for all z. Consequently, if / is a second solution 
of (5.1) then f is rational, entire, and 1 at oo, hence identically 1. Thus the solution 
is unique and m is rational in x. Now the potential q = —[J , fi], so q is also rational. 

Conversely, the equations (5.2) allow construction of rational solutions of the isomon- 
odromy equation. 
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Theorem 5.2. Suppose that B is a constant matrix with trace zero and minimum eigen- 
value —r, and that exp{2niB) = 1. Suppose that the equations (5.2) have a solution. Then 
q = —[J , fi] is rational and is a proper solution of the isomonodromy equation having 
normalized wave function m = l + z~^ f-i_-\-- • ■ + z~'^fr. Both m and q are rational functions 
ofx. 

Proof. In view of the previous proof, all we need to show is that m = l + z~^ fi-\ \-z~'^ fr 

satisfies equations of the form (1.8), (1-9). Let if) = me* = wz^. By assumption, w is 
entire and invertible everywhere. Consider the function 



1-1 dw _i 
ip = — — w — zJ . 
ox 



which is entire. The behavior at 2; = oo is obtained by rewriting the function as 



dx 



m-' = -[J,fi]+0{z-') 



Thus by Liouville's theorem the function is independent of z and is identically q{x) 
~[J J /i]) cind m satisfies (1.8) with this choice of q. Similarly, consider the function 



A{x,z) = 



dtp 



■ d 

z— I wz 



B 



z-^w- 



dw 



w 



-' + B. 



It is regular at 2; = and its behavior at 2; = 00 is seen by rewriting it in terms of m: 



d_ 
dz 



( A 

I me I 



-1 

e m 



dm 
dz 



m 



+ m {z'^ii + xzJ) m-^ = z""/! + 0(z"-^) 



Thus ^ is a polynomial in z with leading term z'^fj,, and m satisfies (1.9). This completes 
the proof. 

We turn now to a closer inspection of the system (5.2), particularly the behavior for 
large x. Note that the coefficients Cj of the expansion of e* have the form 



(5.3) 



1 



Cj{x) = -x^J^+Oix^-""). 



Rescale the problem for large x by setting x^ fj — fj. Then in view of (5.3) the equations 
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(5.2) have the form 



(5.4) 



frVl-r 
[frJ + fr-l]V2-r 



[hirr + . . . + /i J] = -Fi + 0(a;-"), 



The key to the behavior for large x is thus the rd x rd matrix 



(5.6) 



JV2-r 
V2-r 

L 



1 rJ^-'^r 



(r-2) 







1 TS+r-lT^ 

(s+r-1)! '^s 
1 7"s+) — 2t/ 
(s+r-2)!"' 



Since the Vj are constructed from eigenvectors of S, invertibihty of (5.6) is a condition 
on the matrices J, B afone. 

Theorem 5.3. Suppose that e^'^*^ = 1, the least eigenvalue of B is —r, and the matrices J 
and B are such that the matrix (5.6) is invertible. Then the equations (5.2) are solvable for 
all large x and the associated potential q is rational; moreover q{x) = x~^qQ + 0{x~^~'^), 
where the constant matrix qo is similar to B. 

Proof. The preceding discussion estabhshes that invertibihty of (5.6) imphes solvabihty of 
(5.2) for large x. It is clear from (5.6) that the fj are of the form x~^ fj = x^^ fjo+0{x~^~'^) 
with fjo constant. Therefore q = — [ J , /i ] has the form asserted above. Finally, it follows 
by induction that Fj of Theorem 1.1 is 0{x~^). Now B is similar to F^ + xq = qo + 0{x~'^) 
for every x, so B qo. 



Example. As a simple example to illustrate the construction of rational solutions we 
consider the AKNS hierarchy; cf. also [FN]. Here 





"1 


0" 




11 


u 


= (73 = 





1 


Q 


V 
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Prom the recursion relations (1.5) we find 



Fo = as, Fi = q, 



—UV Ux 
—Vx UV 



and the isomonodromy deformation equations (1.6) at order 2 are 

{xu)x + \Uxx - U^V = 0, {xv)x + \Vxx + UV^ = 0. 

Under the symmetry reduction u = ±v these equations reduce to the single equation 

{xu)x + \uxx ± tt^ = 0. 

For the case r = 1 equations (5.3) are simply 

hvi = Q, (l + /iCi)^;i = 0, 

where Ci = xus, and vi is the eigenvector of B with eigenvalue 1. Since we do not know 
B a priori, let us take v = [1, /i]* where h is to be determined. We find 



-1 h-^ 
-h 1 



m = l^-z Vi> g = -[cr3,/i] = 



h-^ 

h 



It is easily verified that det (1 + fiz ^) = 1 for all h. Moreover this result is independent 
of n, the order of the equation in the hierarchy. Thus for all h we get a solution with 



(5.^ 



1 h ^ , 

u = —- — , V — , h 0. 

hx X 



The construction of rational solutions for r = 2 is somewhat more involved. Under the 
symmetry reduction u = v, we have ani{x, z)a~^ = m{x, —z), where 



1 

1 



hence 



For n = 2, 



a, {-lybj 
bj {-lyaj 



X 



Ci = xaa, C2 = yl + 0-3, C3 



-^aa + xl. 
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In the case r = 2, and h = 1 one finds 

/ 3(x^ + l)(x-l) 1\ 
^ - - (, x{x^ + 3) - X J^- 

It is also interesting to consider the symmetry reduction v = eu where e = ±1. Replacing 
t by it in the time evolution equation qt = [-^—Qi F^] , one obtains the nonlinear Schrodinger 
equation 

(5.9) iut = ^Uxx - e\ufu. 

There are solitons only in the case e = — 1, called the self- focussing case. This equation is 
invariant under the scaling u{x,t) —>■ Xu{Xx, X^t) and the associated similarity solutions 
take the form 

1 X 
u{x,t) = ^ = ^ 

where satisfies the ordinary differential equation 

(5.10) z(e0)' + 2e|0p0-0" = O. 

We can construct "rational solutions" of this equation for real x, although they cannot 
be continued into the complex x plane. For example, taking h = i m (5.8) we obtain 

i —i_ 
u — —, V — — = u 

X X 

for real x. Hence we obtain rational solutions of (5.10) for real x in the defocussing case. 

In the self-focussing case, xq+F2 is skew-Hermitian, so always has imaginary eigenvalues. 
Therefore the solutions wz^° have non-trivial monodromy in a neighborhood of the origin, 
and solutions of (1.8), (1.9) for m rational in z do not exist. Thus the theory of the present 
section does not apply. 



6. Backlund transformations 

In this section we adapt a procedure which goes back at least to Moutard [Mt] for 
transforming a pair consisting of a linear differential operator and its eigenfunctions to a 
new such pair. Here the idea is to make a gauge transformation of normalized solutions ip 
associated to a given problem. We look for a gauge transformation of the form 

(6.1) i'ix, z) = G{x, z)ip{x, ^) = [l — z~^a{x)^'i(j{x, z) . 
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We assume that if) is invertible and satisfies (1.1), (1-2), where q satisfies the isomonodromy 
equation. Then ij) satisfies 



(6.2) ^ = [G{Jz + q)G-^ - G^G-^] ^ , 

(6.3) = [GA{z, x)G-^ - zG,G-^] ^ = A{x, z) ^ . 

oz 

In order for 6.2 to take the form dil^/dx — {Jz + g) '0 it is necessary and sufficient that 

da 

(6.4) g = g+[J,a], — = [q + Ja, a]. 



This matrix Riccati equation for a can be hnearized by setting a = bcb~^ with c constant; 
then the equation for a is satisfied provided db/dx = qb + Jbc. In particular, there is a 
solution for each choice of c and each choice of (invertible) initial value b{xo). Note that 
= implies = 0. 

The next problem is to ensure that A{x, •) in (6.3) is a polynomial with leading term 
z'^fi. 

Theorem 6.1. Suppose a satisfies (6.4), anda^ = 0. Suppose thatA(xo, •) is a polynomial 
in z. Then A{x, ■) is a polynomial in z for each x for which A is defined, and q satisfies 
the isomonodromy equation. 

Proof. The assumption — implies that G~^ = 1 + ^"""^a, so 

(6.5) A{x, ^) = [1 - z-'^a{x)]A{x, z) [l + z-'^a{x)] - z-'^a{x) . 
The zero-curvature condition (1.3) is preserved by gauge transformations, so 

(6.6) [^-zJ-q,A]^zJ. 

Equation (6.5) shows that A has highest term z'^/j, and (6.6) shows that if A{x, •) is regular 
at 2; = for one value of x, then this is true for all values of x. To show that (6.6) amounts 
to the isomonodromy equation for q, we must relate A to the coefficients of the expansion 
of rhumT^, where in = V'e* The required relation (1.4) follows from the second part of 
Theorem 1.1: let V'l/ = rn,ye^ be the normalized wave function for the sector fij^. Then 
777-1/ = Grrii, has an asymptotic expansion as 2; — > 00 in ili^ with leading term 1, as required. 
Thus q satisfies the isomonodromy equation. 

We now consider the case of proper q, i.e. q for which the fundamental solutions have 
the form (3.3) ai z — 0. 
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Theorem 6.2. Suppose that q is proper and that xq + is diagonalisable. Then either q 
satisfies an algebraic differential equation of degree n — 1, or there is a nontrivial Bdcklund 
transformation q ^ q of the form (6.4) such that the eigenf unctions transform by (6.1) and 
such that q satisfies the isomonodromy equation. Moreover, q is also proper and xq + Fn 
is diagonalizable. 

Proof. Fix Xq and write matrices with respect to a basis for which Aq[xq) is diagonal. 
Choose a solution of the form (3.3), so w(a;o, 0) = 1. Then detw = 1. Write w{xq,z) = 
Y^z^wjf. Suppose that for some j ^ k, {wi)f~j 7^ 0, and take 

(6.7) djkWijk = 1 ; all other entries of ars{xo) — 0. 

Then a simple calculation shows that the jj-entiy of G{xo, z)w{xo, z) vanishes at 2; = 0, 
and all entries except the jfc-entry are regular at 2; = 0. Therefore 

G{xq, z)w{xo, z) = w{xo, z)z^ , 

with w entire, dettu = 1, C diagonal, C = ejj — ej-k- It follows that 

d 

A{xq,-) =[z—w + wB]w~'^ , B = C + Aq{xq). 

Thus ^(xo, ■) is regular at ^ = 0. We can prolong a by (6.4) and deduce the result from 
the preceding theorem. Clearly q is regular at a^o and hence for all x. 

The obstruction to this construction at x = xq is the vanishing of the off-diagonal entries 
of wi(xo), so [ ^0 ) ] = 0. The equation for w is 

dw 

z—— = A(z, w)w — wB . 
dz 

Substitution of the power series for w into this equation yields Wi = [Aq, wi]+Ai; hence one 
would have [Aq, Ai] = at xq- Thus the obstruction is the vanishing of this commutator 
at every point, which is the algebraic ODE for q of order n — 1: 

[Fn + xq, Fn-i+xJ] =0. 

A consequence of the argument just given is that the matrix a for the transformation in 
Theorem 6.2 can be constructed algebraically, bypassing the differential equation in (6.4), 
by looking at w{x,-) in a basis which diagonalizes Ao{x). Note that xq + Fn is similar to 
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B and hence diagonalizable. Relative to xq + F^, one eigenvalue has been increased by 1 
and a second eigenvalue has been decreased by 1. 

Example. We illustrate the method to construct integer solutions of the Painleve II 
equation by Backlund transformations (cf. also Airault [Ai]); we use the notation in §3. 

The Backlund transformation from the zero solution to the solution for r = 1 coincides 
with the reduced wave function m. In fact, (1.8) is equivalent to the intertwining relation 



mi — zJ 

\ox 



d_ 
dx 



— zJ — q )m 



We work in a basis in which q is diagonal. Note that 



1 
72 



1 -1 
1 1 



Thus in this basis J = —a. 

We look for a gauge transformation of the form (6.1) with a{x) = v{x)T. In order that 
a satisfy (6.4) we must have 

-[(7T,T]=AT, T^ = 
for some scalar A. These equations are satisfied by the choices T = N±: 

N+ = 

The Riccati equation (6.4) is 



'o 


o" 


N_ = 


"o 


l" 


1 












hence 
(6.8) 



v^, = v . 



We use (6.7) to get the initial condition at an arbitrary point xq. Since we are gauging 
from the trivial solution, w = exp{—{xz + z^/3)a}. Therefore, wi{xo) = —xoa, while 
a{xo) = v{xo)N±. In either case, (6.7) implies 



v{xo) = - 



Xo 
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Since xq was arbitrary (other than xq = 0) we see that v{x) = —1/x, and this function 
satisfies (6.8). Thus, as remarked above, v is fully determined by the algebraic condition 
(6.7). 

The corresponding potentials and wave functions are 

q± = [J,vN±] = ±-a3, m± = 1+ —N±. 

X xz 

Let us take the case q = q+ and construct the gauge transformation from r = 1 to r = 2. 
We begin by determining wi. The wave function for is 



m+e 



1 

(xz)-^ 1 



-{xz+z^/3)a 



Expanding the exponential to third order terms in z we find 



m+e = 



where 



Therefore 



" -x' 




"1 


0" 




x-^ 


+ z 







+ 




w 





1 + z 



x' 




-w 


-) 






x^ - 1 


w 




3x 



z 
z-^ 



Wl = 



X 
-w 



in (6.7). Then with k = 2,j = 1 we take 



z-^ 
z 



X 



-X 
-1 



v{x) = 



Sx' 



w 1 — x^ 



This choice indeed satisfies the Riccati equation (6.4), which in this case is 

Vx = \-v . 

X 

The potential and wave function for r = 2 is 



Q = - 



6x^ 2 
+ 



m 



1 — a;^ X 

The gauge transformation G is the first factor of m. 



3x^ 
(l-a:3)2 

1 



1 

{xz)-^ 1 
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7. Scaling, self-similarity, and construction of isomonodromy deformations 

The isomonodromy deformations (1.1) and (1.2) can be derived from the Lax pairs for 
integrable systems by a scahng invariance. This was done for the mKdV equation by 
Flaschka and Newell [FN], working with the 2x2 AKNS system. For the general d x d 
system (1.1) and a given constant matrix fi, the n-th equation in the associated hierarchy 
of equations is 



(7.1) ^^=[J,F,,^+i] = 



d__ 



Here Fq^k is the coefficient of 2; in the formal series of Theorem 1.1. This equation is 
associated to the operator 

Dt,q = I - [z-FqU , 

where denotes the polynomial part of z'^Fq. In fact (7.1) is equivalent to the 

zero-curvature condition 

(7.2) [i?,,,, A,,] = 0, 

which is the compatibility condition for the overdetermined system 

(7.3) £'x,gV' = 0, A,qV' = 0, 
where 



Dx,q = — - zJ -q. 



d_ 
dx 

A solution q = q{x ^t) oi {1 .1) is said to be self-similar if 



(7.4) q{x,t) = qx{x,t) = \-'q{\-'x,\-''t), A>0. 

Note that such a function is uniquely determined by its values at fixed t, say at t = 1/n. 
Define dilations Tx, A > 0, acting on functions of (a;,t, z) by 

(7.5) TU{x,t,z) = f{Xx,X^t,X-^z). 
Let q be self-similar; then 



(7.6) 
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The first of these operator identities is immediate from (7.4). For the second, let mo{x, z) 
be a formal solution of Dx,qm = —zmJ at t = 1/n and extend mo as a function of 
t by m{x,t,z) = ■mo{{nt)~^/'^x, {nt)^^'^z). Then m is invariant under the dilations, so 
the first identity in (7.6) shows that m is a formal solution for each fixed t. It follows 
that Fq = mjimT^ is dilation invariant, which implies the second identity in (7.6). A 
consequence is that the evolution equations (7.1) are invariant under dilation. 

The basic observation is that self-similar solutions of the evolution equation (7.1) cor- 
respond to solutions of (1.6) and conversely. 

Theorem 7.1. If q{x,t) is a self-similar solution of (7.1), then q{x, 1/n) is a solution of 
the isomonodromy equation 



(7.7) 



d 

q, xq + Fn 

ax 



0. 



Conversely, if qo{x) is a solution of (7.7), then q{x,t) — (nt) ^^'^qo{(nt) ^/""x) is a self- 
similar solution of (7.1) with q{x, 1/n) — qo{x). 

Proof. Suppose q is a self-similar solution of (7.1). Differentiating (7.4) with respect to A 
at A = 1, we obtain the Euler equation 

(7.8) .|i+„t| + , = 0. 

At t = 1/n this identity becomes qt — —{xq)x. Substituting this last identity in (7.1), we 
obtain (7.7). 

Conversely suppose that qo{x) satisfies (7.7). Let q be the self-similar extension. By 
the argument leading to the second identity in (7.6), Fq is the invariant extension of F. 
Consequently we only need to check (7.1) or (7.2) when t — 1/n, where it is implied by 
(7.7), (7.8). 

It is worth noting that the scaling argument carries over to the wave functions. Scaling 
invariance Txi/j = ijj is equivalent to the condition that '4>{x,l/n, z) satisfy (1.2). In fact 
scaling invariance implies the Euler equation 

and, at t = 1/n, (7.10) and (7.3) imply (1.1). Conversely, (1.1) and (7.10) for a function 
il)o{x,z) imply that the self-similar extension il){x,t,z) = ipo{{nt)~^^'^x, (nty^'^z) satisfies 
(7.5). 
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Similar considerations apply to the Gel'fand-Dikii flows 

n-2 J 

(7.11) L ^ [[L>'/^]+ , L] , L^Lr,^D^ + Y,u,{x,t)D\ ^ ^ IT ^ 

where A; e N is not divisible by n and [L'^/"']^. denotes the differential part of the pseudo- 
differential fractional power. The coefficients of (L'^/'^)+ are polynomials in the Uj^s and 
their derivatives, with no constant terms. The simplest examples are 

L3 = D"^ + U1D + U0, [Ly\ = D^ + lui, 
These lead to the KdV and Boussinesq equations respectively: 

= [ , L2] = ^Uxxx + ^uux ; 

UltD + Uot = [ [LI^^]+ , L3] {2uox - Uixx)D + {uqxx - \uixxx - \uiUix) . 

The next equation in the KdV hierarchy, 

Ui = Yg-D U -\- -^UUxxx ~l~ 'g'^x'^xx ~l~ "s"^ '^xj 

comes from 

[L^/']+ = D^ + luD^ + fuxD-" + l{5uxx + Su^)D + ^(u^xx + 2uUx). 
A solution L of (7.11) is said to be self- similar if 

(7.12) Uj{x,t) = X^~^Uj{Xx,XH) , A>0, j = 0,1, . . .n - 2 . 

Examples. The simplest cases are the self-similar KdV and Boussinesq solutions which 
satisfy, respectively: 

{xu)x + U + \Uxxx + \uux = ; 
{xUi)x = Uixx - 2wOa; , {xUq)x + Uq = \uixx + fwi^la; " Uqxx ■ 

The next equation in the KdV hierarchy leads to 

(xu)x + U+ YqD^U + ^UUxxx + \uxUxx + ^U^Ux = 0. 
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Note that (7.11) is the compatibihty condition for the system of two scalar equations 

dv 



(7.13) 



Lv{x, t, z) = z^v{x, t, z) 



dt 



The self-similarity condition (7.12) is equivalent to a pair of operator identities analogous 
to (7.6): 

(7.14) [L - z"/] Ta = A"Ta [L - z^I] , (L'=/")+ Tx = X^Tx (L'=/^)+ , 



where 



(7.15) 



Txv{x, t, z) = v{Xx, t, X~^z) 



The self-similarity condition is compatible with the scaling condition 



(7.16) 



v{Xx,X''t,X ^z) ^ v{x,t,z) , A>0, 



on solutions of (7.15). 

We study the scalar equation Lv = z'^v by using the standard procedure to convert it 
to a first-order system. Let 



= (-« 



dv 



' 5a; ' ' ' dx'' 



Then Lv = z^v is equivalent to 
(7.17) 



where 



(7.18) 



1 

1 










Uq Ui U2 



The commutator [ {L^/'^)j^ , L] is a differential operator of order n — 2: 
(7.19) [iL^''')^,L] = Y,w,(^\ . 
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If Lv = z^v then derivatives of (L'^/"')+v with respect to x can be expressed in terms of 
derivatives of D^v of order less than n, i.e. the entries of V'- This leads to an identity 



(7.20) 



where the n x n matrix Gk = Gn,k is a polynomial in z"^ and also in the uj and their 
derivatives of order less than k; (7.20) holds for all such solutions v. 



Examples. The two simplest cases are 



G 



2,3 



73 I 1 



—Ux 2tt 
-2z^u - 2u^ - u, 



G 



3,2 



XX Ux 

2ui 
-3uo + 2uix -ui 
2uixx - uqx -3uo + Ulx 
We differentiate (7.20) with respect to x to find that 

d 






-til 



(7.21) 



dx 



- J^-q, Gk 



where — r is the matrix corresponding to (7.19), i.e. 

0. 



(7.22) 



— r 







_Wo Wi W2 ■■■ 0. 
It follows that the Gel'fand-Dikii equation (7.11) is the same as the matrix equation 

(7.23) qt = r. 

With this preparation we turn to the self-similar solutions of the Gel'fand-Dikii equation. 
The self-similarity condition (7.12) is equivalent to 

(7.24) q{Xx, XH) = X-^d{X)-\{x, t)d{X) , 

where d{X) = diag(l, A, A^, . . . , A^-^). 

If a scalar function v is invariant under the dilations (7.15), the corresponding column 
vector i/j transforms by 

(7.25) ip{x, t, z) = d(X) ip{Xx, X% X'^z) . 
Let P denote the diagonal matrix 

(7.26) " 
The Euler equation equivalent to (7.24) is 

(7.27) 

The following result is analogous to Theorem 7.1 and can be proved in the same way. 



P = ^M(A)]a=i = diag(0, 1, 2, . . . , n - 1) . 
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Theorem 7.2. If the matrix q{x, t) corresponds to a self-similar solution of the Gel'fand- 
Dikii equation (7.11), then q{x,l/k) is a solution of the system of algebraic ordinary dif- 
ferential equations 

(7.28) {xq)a,-\-[P,q]-\-r = 0. 

Conversely, the self-similar extension of a solution qo{x) of (7.28) corresponds to a solution 
of (7.11). 

Examples. The simplest cases are the self- similar KdV and Boussinesq solutions which 
satisfy, respectively: 

{xu)x + U + \uxxx + \uux = ; 
Proposition 7.3. Equation (7.28) is equivalent to the commutator condition 



(7.29) 



— -J^-g,^— - x{Jz + q) - P - Gk 







Proof. This is a straightforward calculation using (7.22) and the identity 
(7.30) [P, J,]-2;— (J,) + J, = 0. 

We show in the next section that (7.29) is the equation for a monodromy-preserving 
flow. 

8. Gel'fand-Dikii equations and isomonodromy 

The equation (7.29) which characterizes self-similar solutions of the Gel'fand-Dikii equa- 
tion (7.11), is the compatibility condition for the system 

(8.1) |^ = ['^^ + ^J^' 

(8.2) z^^[xJ, + xq + P + Gh]i^. 

oz 

The basic result for the forward monodromy problem is the following analogue of results 
of 552, 3. 
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Theorem 8.1. The system (8.2) has a regular singular point at z = 0; the Stokes matrices 
for the irregular singular point at z = oo are preserved under the flow (7.28). 

The fact that the origin is a regular singular point for the system (8.2) is a consequence 
of the fact that Gk is a polynomial in z (in fact a polynomial in z^). To prove that the 
Stokes matrices are constant we must analyze the behavior of solutions of (8.2) as z tends 
to infinity. The argument is similar to that in §2, and needs some preparation. 

Note first that in the trivial case q = there is a fundamental solution of (8.1) having 
the form i/j = A^e^^"^ — d{z)Ae^^'^ where 



U) 



A = 



- 1 


1 


1 - 















J = diag(Q;i,a2 ■ ■ • , "n) 



and the aj are the n-th roots of unity; to fix the choice we let aj = exp{27rij / n) . Moreover 



(8.4) 
(8.5) 



(8.6) 



[ 2; J , m ] + qm . 



It is natural to look for solutions to (8.2) in the form ip — Azme^^"^ , so (8.1) and (8.2) 
become 

dm 
dx 

(8.7) ^ ^^^"^ + xq + G]m - m[[zJ)^ + xzJ] , G = A^^GA^ 

Equation (8.6) has formal solutions as in §1: m = Z^j^o ^ "'/i with /o = 1. 
Example. With n = 2 and with xq fixed, the unique formal solution with diagonal part 
= 1 sX X = xq is 



m{xQ, z) = 1 + {\z u+j^z [uxx + 2u]) 



1 

1 



-1 

1 



+ 0{z-^). 



For later use we note a symmetry property of (8.6). Let 11 be the permutation matrix 



n = Ji = 



1 

1 

1 
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Then simple calculations show 

(8.8) U{zJ)U-^ = {az)J, Uq,U-^ = qaz, Ilf,U-^=faz, UA^U'^ = A^z- 

Note that — A~^ Gk A^ is a rational function of z with leading term z^'J^ as z tends 
to infinity. 

Example. With n — 2 and k — 3, 



Gs = z^J + \zu 



' l" 




"0 


1" 


-1 


4''^ a; 


1 






1 1 
-1 -1 



The next step in our analysis of the behavior of solutions for large z involves under- 
standing Gk- 

Lemma 8.2. Ifm is a formal solution of (8.6), then the product m{zJ)^ differs from 
Gk only in terms of negative degree in z. 

Proof. Equations (7.21) and (8.6) imply that 

d 



dx 



— zJ , rh ^Gkrh 



We write 
(8.9) 



m 



3^0 



771 ^rrh = 0{z ^) 



and take gj = hj = for j < 0, so that the equation above is equivalent to 
(8.10) [J , gj+i ] = {gj),: - hj-k, all j. 

It follows from the relations (8.10) and the identity go = that gj is diagonal and constant 
for i < k. Note that this fact does not depend on (8.2); it is true for arbitrary smooth q. 
As in §1, the fact that the diagonal part of rri~^Gk'fn is independent of the choice of formal 
solution rh implies that the coefficients of the diagonal part are polynomials in the entries 
of q and their derivatives, and that all but the top order term (zJ)^ are polynomials with 
constant term zero. The only such polynomial which is constant for arbitrary q is the zero 
polynomial. Therefore gj — for < j < k and the assertion is proved. 

As in §2 we may use a partial sum of a formal solution to regularize our problem at 
infinity. Let / = Yl'i=o fj ^® ^ partial sum of m; we have 



f-'Gkf^z''j'' + 0{z-'). 
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We look for a solutfon of (8.2) in the form ij) = f ip. Note that 

(8.11) .|^-pV = a..^(a;V) 

so the equation for -ij) is 

(8.12) z ^ = [-zf-'f, + xf-\zJ + q)f + r'Gkf] V' 

= [z'^j'' + xzJ + r{x, z)] ip , r{x, z) = 0{z~^) . 

From this point the analysis is essentially the same as in §2. Equation (8.2) has unique 
solutions of the form 



(8.13) V^^ = A^m^e*, <^{x, z) ^ ^z'' j'' + xzJ , 



where has an asymptotic expansion with leading coefficient 1, valid in a sector fl^. As 
in §2 these functions also satisfy (8.1) and therefore the Stokes matrices which relate the 
i/^i, are constant. This completes our sketch of the proof of Theorem 8.1. 

Remark. The symmetry property (8.8) and the uniqueness of the solutions (8.13) allow 
us to conclude that these solutions have the n-fold symmetry which can be stated roughly 
as 

m(x, az) = Ilm(x, z) 
To state it more precisely we need to investigate the set 

T,= {z : Re{zaj)'' = Re{zai)'' , some ^ af}. 

Let p be the smallest positive integer such that a'^'' = 1, i.e. such that n divides pk. Then 

'E = {z:z''eiR}ifp = 2; S = : z'^p'' e R+} if p > 2. 

We define regions Qi, as in §2, bounded by various of the rays of the set E. It follows from 
the form of S that fl^+i is obtained by rotation of through an angle tt /pk, {ir/k if 
p = 2) and the precise form of the symmetry is 



(8.14) mt,j^s{xi(y.z) = Tlmv{x,z)'n. | 



s = 2kp/n if p > 2 
s = 2k/n ii p — 2. 
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The Stokes matrices satisfy the corresponding symmetry 

(8.15) S^+s = US^U-\ 
as weU as the standard constraints 

(8.16) diag(S'j,) = 1, e^S^e~'^ is bounded as 2; — > oo, z & ilj^, z near Sj^, 

where T,^ is any ray in 0,^, fl fij^+i; it will be convenient to choose Ej^ to bisect this sector. 

The forward monodromy problem at 2; = is the same as in §3, but we must take the 
symmetry into account. For convenience we consider only the generic case, when there as 
a traceless constant matrix A and a solution of ip of (8.1), (8.2) with the property that 
'ip{x, z)z~^ is an entire function of z, invertible for each z. If so, then there are traceless 
constant matrices Bi, such that '^i,z~^'' are entire, and 

(8.17) B^+i=S-^B^S^; exp{27T iB^) ^ S^S^+i . . . S,-i, B^+s = UB^U-\ 

Theorem 8.3. In the generic case, if two solutions of the isomonodromy equation (7.29) 
have a common domain and the same monodromy data {Si,, B^}, then they are identical. 

Proof. Suppose that q and q' are two solutions of (7.29) on a common domain, having 
the same monodromy data. Let i/'i/ and tj)'^ be the corresponding normalized solutions of 
(8.1), (8.2). Let m = 'i/ji,e~'^ for 1^1 > 1, 2 between Ti^-i and Ej^ and m = 'i/jiyZ~^''' e~'^ on 
{1^1 < 1}. Define m' analogously. Each of m and m' is a solution of a matrix Ricmann- 
Hilbert factorization problem; a vector version of this problem is described in more detail 
below. It follows that m~^m' is continuous and piecewise holomorphic, hence entire. From 
the asymptotics of (8.13) we know that A~^m ~ 1 at infinity, and the same for m', so 
Q = m~^m' is a polynomial in z and 

d{z)-^Qd{z) = l + 0{z-^). 

This last fact implies that Q — 1 is strictly lower triangular. Therefore m' = Qm has the 
same first row as m. But m and m' are each determined by the first row M; in fact the 
j-th row is D^+^(Me^) e"*. Thus ip = ip' and q = [DiP_J^iP]ip-^ = q' . □ 

In order to formulate the inverse problem, in which matrices {Si,,B^} are given, we 
must give a precise description of the vector Riemann-Hilbert problem alluded to above. 
Let M be the first row of 'ipi^e~'^ for \z\ > 1, z between Ej,_i and E^. On {\z\ < 1} 
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let M be the first row of 'ipi,z~^"e~^ = h-zfriy; this is independent of v. Because of the 
asymptotics of rrii, and the relations given by the Stokes matrices, one can see that the 
row vector function M has the properties 

(8.18) M has limit (1, 1, . . . , 1) as 2; ^ 00; 

(8.19) M is bounded and holomorphic where defined; 

(8.20) M is continuous up to the boundary from each component; 

Moreover the boundary values M^^^ the limit on the boundary of the region outside the 
unit circle between the rays and Sj^, and Mr, the limit on the unit circle from the 

disc, satisfy 

(8.21) M^+i = M^e^S^e-^; = Mre^z^-^e'^, 
while the row vector M itself has the symmetry 

(8.22) M{x,az) = M{x,z)U-^, \z\>l, z ^[Je^. 

Conversely, suppose matrices {S,^, By} are given. Generically the Riemann-Hilbert prob- 
lem (8. 18)- (8.22) has a unique solution, say for x in some domain. 

Theorem 8.4. Suppose that S^, and B^, are constant matrices which satisfy the conditions 
(8.15)-(8.17), and trBi, = 0. Suppose the Riemann-Hilbert problem (8.18)-(8.22) has a 
unique solution M{x, ■), forx in some domain. Then there is a unique n-th order operator 

(ajJu 

such that the row vector v = Me* satisfies the equation Lv = z^v. The system correspond- 
ing to L satisfies the isomonodromy equation (7.29). 

Proof. The vector function M has an asymptotic expansion in powers of 2;"-'^ as 2; — > 00, 

CXD 

(8.23) M{x,z)r^Y^ z ^aj{x) 

and the symmetry (8.22) implies that the coefficients have the form 



(8.24) 
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where the 6j's are scalars. The expansion can be differentiated term by term with respect 
to X. 

Next we construct the sequence of vectors 

(8.25) M(^) = D^(Me*)e-* = (L»^v)e-*, j > 0. 
Note that 

(8.26) M^^) = (1,1,..., l){zjy + 0{z^-^). 

In particular M^"') — z'^M is 0{z'^~^). We use the symmetries (8.24) and their analogues 
for the M^^^^s to conclude that there is a unique scalar function Un-i{x) such that 

(8.27) M(") - + «n-iM(^-i) ^ 0(^^-2). 
Continuing, we find unique functions Un-2, ■ ■ - Uq such that 

n-1 

(8.28) M(") - ^"M + UjM^^^ = 0{z-^). 

Now the M(^)'s and z'^M are solutions of the Riemann-Hilbert problem (8.19)-(8.22). 
Therefore the left side of (8.28) is a solution, with limit as 2; tends to 00. By our 
uniqueness assumption for solutions of (8.18)-(8.22) the left side of (8.28) must vanish. It 
follows that the row vector v = Me* satisfies the equation 

n-1 

(8.29) Lv =[D'' + J2 UjD^]v = z'^v. 

j=0 

Our next task is to show that Un-i = 0, so that L has the desired form. 

We may reduce to a first-order system as before and conclude that v is the first row 
of a solution to a system of the form (8.1). The j-th. row of the matrix ip is precisely 
jVf(i-i)e*. Note that the trace of the matrix q is —Un-i- It follows that the determinant 
A = det(V') = det(V'e~*) satisfies 

(8.30) A^ = -un-i A . 

On the other hand, A is a solution of a scalar Riemann-Hilbert problem with trivial 
multiplicative jumps. Therefore A is an entire function. The asymptotic behavior is the 
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same as that of detA^, so A is a polynomial with leading coefficient detA^. Since this 
coefficient does not depend on x, (8.30) implies Un-i = 0. 

Similarly, zM^ - xM^^^ ^M[{zJf + xzJ] - M^^'^ is a solution of (8.19)-(8.22) which is 
0{z^~^) and we may use the symmetries and (8.26) to deduce that 

zM^ - xM^^^ + M[{zjf + xzJ] = M^'^) + au-iM^^-^^ + 0{z^-^). 

As before, we conclude that v satisfies an equation of the form zVz = xv^ + L^w, where 
L* is a differential operator in the x-variable, of order k. This equation is equivalent to 
the first-order system 

= a;(J^ + q)i; + P^j + G^j, 
where the matrix G{x,z) is characterized by the condition analogous to (7.20): 

(8.31) ( L*v, DL*v, . . . , D^'-^L^v Y = Gij;. 

Equation (7.29) is the compatibility condition for (8.31) and (8.31), provided we may 
replace L* by {L''/'') + . Wc may make this replacement if and only if L^v — (Zy'^/"^)_|_f , 
and this is true if the matrices G and Gk have the same first row. 

We write ip = K^me'^ and note that m has an expansion at oo with leading term 1. 
Moreover, m satisfies 

xrux = [xzJ , m] + xqm; zrriz = [xzJ , m] + [xq + G]m — m{zJ)^, 

where G = A~^GAg. It follows from this and Lemma 8.2 that 

(8.32) G = m{zJ)^m-^ + xm^ - zm^ = m{zjfm-^ + O^z'^) ^ Gk + 0{z-^). 
Therefore 

d{z)-'[G - Gk]d{z) = A[G - Gk]A-' = 0{z-') 

which implies that G — Gk is strictly lower triangular. In particular G and Gk have the 
same first row. This in turn implies that L^v — (L'^/")_|_f , so we may replace by 
(L'=/") + . □ 

9. Isomonodromy deformations and string equations 

In this section we consider a second classes of (systems) of ordinary differential equations 
associated to the Gel'fand-Dikii hierarchy (7.11) and show that each is an isomonodromy 
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deformation for a first order system. This class was obtained by M. Douglas [Do] in 
connection with theories of quantum gravity. The Gel'fand-Dikii flow itself is replaced by 
the equation 

(9.1) [(L^/")+,L] = n/ 

where / is the identity operator. For convenience of notation we consider h= 1. The case 
n = 2, k = 3 and the case n = 3, k = 2 each lead to a form of the PI equation; for n = 2, 
k = 3 the equation is 

for n = 3, A; = 2 the constants are different. The case n = 2, A; = 5 is 

j-gZ) U -\- "^UUxxx ~l~ "^^x^xx ~l~ "g"^ ^x ~l~ 1 — 0* 

Remark. The development in this section applies equally to the case when the operator 
(L'^/'^)-!- is replaced by a finite linear combination ak{L^^^)+ with arbitrary constants 
afe. The phase function # which is considered below must be changed accordingly. 

Note that (9.1) can be written as a comutator equation 

(9.2) [L-^"/, ^ -nz"-i(L'=/")+] =0, 

which is the compatibility condition for the system of two scalar equations 

(9.3) Lv = z''v, nv,^nz''{L''/'')+v. 
Reduction to a first order system as in §8 gives a matrix version of (9.2): 

(9.4) ^ - (j^ + (j) ^ z ^ - nz^'Gk =0. 



d d 
^ - (^z + 9) , ^ ^ - nz^'Gk 



In turn, (9.4) is the compatibility condition for the system 

(9.5) ^ = [j^ + ,]V,, 

(9.6) z^^nz^Gki^. 

We sketch a proof of the following analogue of Theorem 8.1. 
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Theorem 9.1. The system (9.6) has a regular point at z = 0; the Stokes matrices for the 
irregular singular point at z = oo are preserved under the flow (9.4)- 

Since Gk is a polynomial in z^ it follows that the origin is a regular point for (9.6). 
The second assertion is a consequence of Lemma 9.3 below. We look for exact solutions 
of (9.6) having the form ij; = A^/i/', where / = ^^^q^^ z~^mj. Then as before there are 
solutions t/ji, of this form which have the appropriate asymptotic expansions in sectors 
and which are related by constant Stokes matrices S^. 

Thus we look for solutions of (9.5), (9.6) having the form i/j = A^me*^ for some diagonal 
matrix- valued function $. The equations become 

(9.7) DxTn = rrix — zJm — qzm = —m^x 

(9.8) Dziri = zrriz - nz^Gkm + A~^PAzm = -zm^z ■ 

Lemma 9.2. Suppose m is a formal solution of (9.7). Then the diagonal part of the 
coefficient gj of z~^ in the formal series rh~^rzrh has the form CjJ~^, where Cj is a 
scalar function which is a polynomial without constant term in the entries of q and their 
derivatives. Moreover, Cj = if j is divisible by n. 

The proof is postponed. 

Lemma 9.3. Suppose (9.1) is satisfied. Then there is a unique diagonal matrix-valued 
function 

n-l 

(9.9) ^x, z) = ^{zJr+'' + J2 H^jy +xzJ-\{n- 1) log^l , 

i=i 

where the bj's are scalar constants, and a unique formal power series in z~^ , m = 
^°^Q 2~-'mj(x) such that mo = 1 and m is a formal solution of the equations (9.7), 
(9.8). 

Proof. We begin by refining Lemma 8.2. Let m be a formal solution of (9.5) and use the 
notation of (8.9). 

We already know from the proof of Lemma 8.2 that gi = ■ ■ ■ = g^ = It follows from 
Lemma 9.2 and the recursion relations (8.10) that the diagonal part of g^j^j has derivative 
CjJ~^ for j > 0. The diagonal part of g^+j itself is a polynomial without constant term in 
the entries of q and their derivatives, so we can write 
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where bj is scalar polynomial in the entries of q and their derivatives, without constant 
term. 

Equation (9.1), which we now invoke, implies that wq = —1 and Wj = 0, j > 0. 
Therefore hj = for j < n — 1 and 

We use the recursion relations (8.10) again and conclude that Qk+j is diagonal and is 
constant for < j < n — 1. Thus the scalar bj is constant, i.e. this expression in q and its 
derivatives is invariant under the x-flow (9.1). 
At the next step (8.10) implies 

gk+n-i = ^{xJ + biJ), bi constant, scalar; 

{9k+n)x — [J J gk+n+1 ] + ^n- 

By Lemma 9.2, hn is off-diagonal and we find that the diagonal part of gk+n is constant. 
This fact is independent of (9.1) and holds for arbitrary q, which again allows us to conclude 
that the diagonal part of gk+n vanishes in all cases. 

Note that the diagonal part of A~^PAg = A~^PA is |(n — 1) 1. The preceding argument 
shows that 

(9.10) 

m-^[nz''Gk - A-^PA^]m 

n-l 

= nz^^+^J^ + ^ jbjizjy +xzJ - ^{n- 1)1 + E 

3 = 1 

p. n—l 

= [^TT^^"^' + E ^j^'^jy +^^J-\{ri- 1) log^l] + E 
= z^^[^{x,z)]+E 

where the bj are scalar constants and $ is defined by (9.9). The remainder term E is 
a formal series with terms of non-positive degree in whose term of degree in z is 
off-diagonal. 

Finally we look for a solution to (9.8) in the form m = ph where p is the sum of the 
terms of degree > —n — k in the formal solution rh. An analogue of (9.10) holds with p in 
place of m, so (9.8) is equivalent to 

zh,^[^,,h]-Fh, F = 0{z-^) 
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where the off-diagonal part of F is 0(2 ^). The proof of Theorem 2.1 in [CL, Ch. 5] shows 
that this equation has a formal solution h = l + 0{z~^). 

We return to unfinished business. 

Proof of Lemma 9.2. The assertion is independent of the choice of formal solution, so 
we may take rh to be the unique formal solution normalized so m{xQ.,z) = 1. Because 
of uniqueness and the symmetry properties (8.8), we can conclude that nm(a;, 2;) n"-*^ = 
m(a;, az). With this and another use of (8.8) we find that 

nm(a;, z)~^rzrh{x.i z) = rh{x, az)~^razTri{x, az). 

Therefore the coefficients in the expansion m~^rzm = ^ z~^hj satisfy 11 hj 11"^ = a~^hj. 
This in turn implies that 

U{J^hj)U-^ = {UJ^U-^){UhjU-^) = {a^J^){a-^hj) = Phj. 

It follows that the diagonal part of J^hj commutes with 11 and therefore is a scalar Cj. 
We noted earlier that since the diagonal part of hj is independent of the choice of m, its 
entries are polynomials in the entries of q and their derivatives. 

When j is divisible by n, the diagonal part of hj has trace ncj, so to complete the proof 
it is enough to show that the trace is zero. But 

tr{m~^rzm) = t'r{rz) = tr{A~^rAz) = tr(r) = 0. 

Examples. 1. In view of the previous calculations for the case n = 2, A; = 3, we find that 
in general 

Z"'^ {Uxxx + QuUx) -2Uxx-Qu'^ 

Equation (9.1) in this case is u^xx + ^uux = —4, so Uxx + ^u^ = —4:{x + b) for some constant 
b. Thus the general formula reduces to 

m-^{2zGz)m = 2z^J^ + {x + h)J - ^z'^ 

It follows that for this case 

^{x,z) = lizJf + {x + h)zJ - ilog^l. 



rh~^{2zG^)rh = 2z^J^ + \ 



1 

1 
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2. For the case n = 3, k = 2, in general 

Equation (9.1) is equivalent to: u\xx — 2tto = 2^2 and |wia;x + \u\ — uqx = hi-\-x for some 
constants 62, hi. Therefore 

$(a;,^) = K^J)^ + 62(2;J)^ + (pi + x)zJ -\ogzl. 

Remark. The set of solutions of (9.1) is invariant under translation, so one can always 
get rid of the constant hi in $ by translating. 

Remark. As wc noted earlier, the origin is a regular point for the equation (9.6), which 
is satisfied hj if) = A^me*. Therefore me* has trivial monodromy. It follows that the 
product of the Stokes matrices is (— 1)"'~^1: 

(9.11) SiS2---So = {-lT-^l. 

The following is the analogue of Theorem 8.3, and is proved in the same way. 

Theorem 9.4. If two solutions of the isomonodromy equation (9.4) have a common do- 
main and the same monodromy data {Si^}, then they are identical. 

The direct problem is similar to, but somewhat simpler than, the direct problem treated 
in §8. The set S has the form 

(9.12) E = : e zR} if p = 2; S = : ^2pn+2pfe ^ if p > 2, 

where again p is the smallest integer such that n divides pk. The symmetry condition (8.8) 
and its consequences carry over to the present situation, with s = 2p{n + k)/n if p > 2 or 
s = 2{n + k)/n ii p = 2.. 

Again we associate a Riemann-Hilbert problem. Let '^i. = Agm^e^ with nii, — l+0(^~^) 
in Qi, and let M be the first row of mj/ for \z\ > 1, z between T,i,-i and 2,^; on the unit 
disc let M be the first row of ^("-^^/^V^oe"*. Then M satisfies (8. 18)- (8.20), and (8.22), 
and its boundary values satisfy 

(9.13) M^+i = M^e^S^e-^, \z\ > 1; 

(9.14) M^+i = z-^"-i)/2^r e^SoSi ■ ■ ■ S^e'^ , \z\ = 1. 

For the inverse problem, suppose that the Si, and $ are given. Generically the Riemann- 
Hilbert problem (8.18)-(8.20), (9.13), (9.14) has a unique solution M{x,-) for x in some 
domain. 
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Theorem 9.5. Suppose that the constant matrices satisfy (8.15), (8.16), and (9.11) 
and that the matrix ^{x, z) has the form (9.9). Suppose that the Riemann-Hilhert problem 
(8.18)-(8.20), (9.13), (9.14) has a unique solution for x in some domain. Then there 
is a unique n-th order operator L = D'^ + Yl^=o '"i-^"' such that v = Me* satisfies the 
equations Lv = z'^v, zvg = nz'^{L^/'^)^v . In particular, the isomonodromy equation (9.1) 
is satisfied. 

Proof. Arguing exactly as in the proof of Theorem 8.4, we examine the expressions 

- ^"M, zM^ + M{z^^) - nz'^M^^'^ 

and show that v = Me* satisfies an equation of the correct form Lv = z'^v as well as 
a second equation of the form zv^ = nz^L'^v. Here L# is an operator of order k. The 
corresponding first order systems are 

'ipx = [Jz + q]^, zipz = nz'^Gip, 

with G determined by (8.31). As before we write ip = A^me*, so m = 1 + 0{z~^) and 

zmz = [nz'^G - A~^PA^] m - 'm{z^z) = [nz'^G - A~^PA] m - m(^$^). 

Because of the form of $ and the asymptotic behavior of m, this equation implies 

G = m{zJfm-^ + 0{z-^) = + 0{z-^). 

As in §8 we conclude from this fact that G and Gk have the same first row and therefore 
that we may replace by (L'^/"')+. 

References 

[ASl] M. J. Ablowitz and H. Segur, Exact linearization of a Painleve transcendent, 
Phys. Rev. Lett. 38 (1977), 1103-1106. 

[AS2] M. J. Ablowitz and H. Segur, "Solitons and the Inverse Scattering Transform", 
SI AM Studies in Applied Mathematics, Philadelphia 1981. 

[Ai] H. Airault, Rational Solutions of Painleve Equations, Studies in Applied Mathe- 
matics 61 (1979), 31-53. 

[BC] R. Beals and R. R. Coifman, Scattering and inverse scattering for first order 
systems. Comm. Pure Appl. Math. 87 (1984), 39-90. 



49 



[CL] E. A. Coddington and N. Levinson, "Theory of Ordinary Differential Equations", 
McGraw-Hill, New York 1955. 

[Do] M. Douglas, Strings in less than one dimension and the generahzed KdV hierar- 
chies, Phys. Lett. 238B (1990), 176-180. 

[FN] H. Flaschka and A.C. Newell, Monodromy and Spectrum Preserving Deforma- 
tions. I, Comm. in Math. Phys. 76 (1980), 65-116. 

[FZ] A. S. Fokas and X. Zhou, Integrabihty of Painleve transcendents, to appear. 

[Ga] E. Gambler, Sur les equations differentielles du second ordrc et du premier degre 
dont I'integrale generale est a points critiques fixes. Acta Math. 33 (1910), 1-55. 

[IN] A. R. Its and V. Y. Novokshenov, " The Isomonodromic Deformation Method 
in the Theory of Painleve Equations," Lecture Notes in Mathematics no. 1191, Springer 
Verlag, Heidelberg 1986. 

[JMU] M. Jimbo, T. Miwa, and K. Ueno, Monodromy preserving deformation ofhnear 
ordinary differential equations with rational coefficients, Physica D 2 (1981), 306-352. 

[JMl] M. Jimbo and T. Miwa, Monodromy preserving deformation of linear ordinary 
differential equations with rational coefficients, fl, Physica D 2 (1981), 407-448. 

[JM2] M. Jimbo and T. Miwa, Monodromy preserving deformation of linear ordinary 
differential equations with rational coefficients. Iff, Physica D 4 (1983), 26-46. 

[Ma] B. Malgrange, "Equations Differentielles a Coefficients Polynomiaux," Birkhauser, 
Boston 1991. 

[Mi] T. Miwa, Painleve property of monodromy preserving equations and the analyt- 
icity of the r function, Publ. R.I.M.S. Kyoto University 17 (1981), 703-721. 

[Mol] G. Moore, Geometry of the string equations. Comm. Math. Phys. 133 (1990), 
261-304. 

[Mo2] G. Moore, Matrix models of 2D gravity and isomonodromic deformation. 
Progress Theor. Phys., Suppl. no. 102 (1990), 255-285. 

[Mt] T. F. Moutard, Note sur les equations differentielles lineaires du second ordre, 



50 



C. R. Acad. Sci. Paris 80 (1876), 729. 

[Sa] D. H. Sattinger, Hamiltonian hierarchies on semisimple Lie algebras, Studies in 
Applied Math. 72 (1985), 65-86. 

[Zu] V. Zurkowski, "Scattering for first order linear systems on the line and Backlund 
transformations," Ph.D. Thesis, University of Minnesota - Minneapolis 1987. 



